Reverse Rotations in the Circularly-Driven Motion of a Rigid Body 
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We study the dynamical response of a circularly-driven rigid body, focusing on the description 
of intrinsic rotational behavior (reverse rotations). The model system we address is integrable but 
nontrivial, allowing for qualitative and quantitative analysis. A scale free expression defining the 
separation between possible spinning regimes is obtained. 
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Nontrivial effects may arise in the realm of integrable 
classical mechanics. The most recent manifestation of 
this fact is the refreshed interest in the dynamics of slid- 
ing bodies In particular, the work by Farkas ei 0, 
an investigation on the frictional coupling between trans- 
lational and rotational motions, has deserved attention, 
being followed by a dozen of papers on related topics 
0. Most of these previous work address the influence 
of friction in the free (not driven) dynamics of disks. In 
the present paper we deal, in a sense, with a comple- 
mentary problem, namely, that of frictionless dynamics 
of a forced rigid body (RB). More specifically, we study 
the dynamical response of RB's submitted to rotational 
driving forces (in a way that will become clear below), 
focusing on the appearance of reverse rotations. Gener- 
ally speaking, a reverse rotation occurs when a system, 
or part of it, is forced to rotate counterclockwise and its 
intrinsic angular degree of freedom develops a clockwise 
motion, or vice- versa. This terminology has been recently 
used in the literature to designate unexpected rotational 
behavior of a cylinder inside a rotating drum filled with 
a viscous fluid Related, though more intricate phe- 
nomena have been reported in the chaotic response of 
a parametrically excited pendulum Q, and in tests of 
printing machinery of journals (|]. Therefore, reverse ro- 
tations are quite a general behavior in diverse physical 
systems. In what follows we propose a simple mechani- 
cal model to investigate this effect in an analytical way. 

Our model system is schematically shown in Fig. [IJa). 
It consists of a uniform disk of mass m and radius R 
resting on a horizontal frictionless surface. We consider a 
disk only for convenience, the following reasoning is valid 
for an arbitrary RB (this point shall be discussed at the 
end of the manuscript). The system is then submitted 
to an external horizontal force F, provided by a driving 
apparatus, through a thin rod attached to a fixed point 
(P) on the disk, around which the whole body can rotate 
freely. The driving apparatus takes the disk from rest and 
makes the point P follow a uniform circular trajectory of 
radius d around a fixed origin (O) with angular frequency 
Lo (see Fig. [Ijb)). For definitencss we assume the rotation 
to be counterclockwise, and without loss of generality we 
use a coordinate system such that the point P lies in the 
positive a;-axis at t = 0. For later times we denote the 
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frictionless surface 



FIG. 1: (a) Circularly-driven motion of a disk on a frictionless 
horizontal surface, (b) The point P, located by the vector d, 
describes a circular path around the origin O. The vector r 
gives the position of the center of mass {CM) and 1 is the 
vector that connects CM to P. 



position vector of P by d and the vector locating the 
center of mass {CM) by r. Since the disk is assumed to 
be perfectly rigid, P is always a distance I apart from 
CM . The relative position of these two points is given 
by the vector 1, as shown in Fig. [TJ^b). Finally, the angle 
between the x-axis and the line connecting CM and P is 
denoted by (j). The variables r and completely specify 
the position of the disk. 

The main goal of this work is to answer the following 
question: Which range of initial angles and param- 
eters m, uj, R, /, and d, lead to a clockwise rotation of 
the disk around P, if this is possible at all? The well- 
posedness of this question shall become evident soon. 
Since the point P itself is being forced in the counter- 
clockwise direction, such a dynamical response charac- 
terizes a reverse rotation (decreasing (/)(t) on average). 

According to our definitions, we have d — dcos{ojt) x-l- 
dsin{ujt) y, where d can take any fixed value in the in- 
terval (0, oo). The vector 1 linking CM to P is given by 
1 = / cos X + / sin 4> y, with < I < R. The position of 
CM is denoted byr = a;x-|-yy. These three vectors 
must satisfy r+1 = d, which provides two holonomic con- 
straints: X — dco8{ujt) — lcos(l), and y ~ dsm{Lut)~lsm(l). 

One initially cope with the problem without being con- 
cerned with the question of how the disk came from rest 
to motion. This point shall be addressed later. Sup- 
pose that Fc, unknown a priori, is the constraint force 
that keeps the circular trajectory of P. This force is 
assumed to be provided by a robust apparatus in the 
sense that the circular path is not affected by the iner- 
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tia of the disk. The equations of motion for the CM 
degrees of freedom read F^^x = "^i, Fc^y — my. The 
Newton's Second Law for the angular coordinate is given 
by Tc = 1 X Fc = IcM(i>2, i.e., 



ml {cos (j) y ■ 



IcM<t> 



(1) 



where Icm is the inertia moment of the disk with respect 
to the center of mass Q . By using the holonomic con- 
straints we find x = I sin</) (f>-\-l cos 4> (jy^ — duj'^ cos{Lot) and 
y = — Zcos(/) -1- /sin0 0^ — duj'^ sin(wi), one can, there- 
fore, decouple the angular equation of motion, which 
becomes: ml[dLu'^ {sin (j) cos{ujt) — cos (f> sin{ujt)) — — 
Icm4>, or 



zldio^ 
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sin(0 — wt) — , 



(2) 



where Ip — Icm + mP . As expected, the uniform plane 
rotation generates an effective gravity. We now proceed 
to the following change of variables: 9 — (p—tut+Tr, imply- 
ing = (p — Lu and 9 — (j). This change transforms Eq. @ 
into the simple pendulum equation 6 + ^j^ui^sin9 = 0. 
Thus, we see that the time dependence of 0, the physi- 
cally relevant variable, is given by a combination of pen- 
dular and uniform motions (plus a constant additional 
factor). 



pendulum 



lu7n{t) -\-LOt -7T 



(3) 



Note that solution ([3]) automatically yields x{t), y{t), and 
Fc{t). It also implies that we have a hidden constant of 
the motion corresponding to the mechanical energy of the 
auxiliary pendulum (we shall use this terminology to refer 
to the pendular term in the solution ^) £ — \lpu + 
2mdlLu^ sin^ (6/2), where the potential energy is set to 
zero in its lower position. In terms of initial conditions 
of the original variable we have 



w)^ -I- 2mdluj'^ cos^ {410/ 2) 



(4) 



By inspecting solution Q we note that reverse rotations 
of the disk are possible only when Openduium describes a 
motion of rotation instead of libration, that is, the me- 
chanical energy of the auxiliary pendulum must satisfy 
£ > 2mdluj^ in negative pendular cycles. Relation ^ 
together with the equality £ = 2mdluP' define the curves 
that separate libration and rotation of variable 9 in the 
space of initial conditions (0Oj 4'o)- One finds that initial 
conditions in between the curves 
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, , ^ . mdl . , , 
l±2xl — sin{(j>Q/2) 



lead to libration of 0, and thus, to normal rotation of 
the disk. Of course, not all points that lie outside this 
region lead to reverse behavior. The pendular rotation 
must be clockwise with an angular frequency larger than 



w, in order to make (f){t) a decreasing function of t, on 
average. The condition for reverse rotation then reads 



T 



for clockwise pendular 



cycles of period T, where K denotes the complete ellip- 
tic function of first kind. In the regions where T equals 
2t:/lu, the constant £ satisfies the transcendental equa- 
tion 
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(5) 



Using this prescription and Eq. ^ we manage to select 
the initial conditions ((/)o, ipo) that lead to reverse and 
normal rotations. These regimes are separated by the 
curve 



Amdlu 
Ip 



■cos2(0o/2) 



(6) 



where £ stands for the solution of Eq. The root 

with a plus sign was discarded because it is related to 
positive pendular rotations. We name the above defined 
curve a synchronization line because initial conditions 
on it develop neither normal nor reverse rotation, for the 
angular frequency of 9 coincides with uj and, therefore, 
the time evolution of </> averages to a constant value. 

At this point we address the problem with static 
initial conditions. This is the most interesting sit- 
uation, since it is expected that high enough clock- 
wise initial velocities 00 trivially lead to reverse rota- 
tions. We note, however, that i, y, and 4> never van- 
ish simultaneously, and, thus, the kinetic energy of the 
disk K. = mx^/2 + my^/2 + Icm'P'^/^ = md'^uj^/2 + 
{IpLo/2 + mdlu) — £ jio) (f)—Ip(tP'/2+Ip(f)'^/2u} , is nonzero 
for all times. It is clear that the force Fc alone is not 
compatible with static initial conditions. In order to en- 
compass these conditions, we assume that an impulsive 
force acts on the disk, taking it from rest to motion in a 
time scale much shorter than any other in the problem, 
e. g. 27r/w. This is a realistic assumption when we have 
a table top engine driving a light body. More explicitly, 
we suppose that the force can be split into two parts 



Fo for t e [0-,0^ 
Fc for t > 0+ 



(7) 



where Fq denotes the impulsive force that acts during an 
arbitrarily small time interval centered at t = 0. In the 
limiting case we have a delta function, which we initially 
write in generic form as Fq = 6{t){a x+/3 y), where a and 
/3 have dimension of momentum. These constants are to 
be determined by the motion we know the force Fg causes 
to the disk. More specifically, we know that, immediately 
after its application, to fit the Fc prescription, the point 
P must acquire a velocity 



vp(0+) = Ludy 



(8) 
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The equations of motion in the infinitesimal interval 
[0~,0+] are mx ~ aS{t), my = f}S{t), and lcoa(j)f3S{t) — 
lsin(j)aS{t) — Icm4>- Integrating one obtains the veloci- 
ties soon after the application of Fq 



i(0+) = ^ , ^(0+) = I , <^(0+) = -i- 

m TO IcM 



(/3cos 00— a sin 0o) 
(9) 



where, here, the static initial conditions were employed: 
i(O-) = 0, j/(0-) = 0, and 0(0") = 0. We also used 
(/)(0^) = 0(0+) = 00 (since the impulsive force causes no 
discontinuity in coordinates). The velocity of P at any 
time is given by vp = r + 1 = (i — /0sin0)x + (y + 
14) cos (f) y. Therefore, the initial velocity as a function of 
a and /? is 



vp(0+) = 



a 
m 



m 



IcM 



sin 00 (/3 cos (po — a sin 0o ) 



■ cos 00 (/3 cos 00 — a sin 0o ) 



By applying the consistency condition ([5]) we obtain a 
pair of equations, involving a and /3, whose solution is 



uidm I 



272 



sin 00 cos 00 , /3 = mduj 1 



ml 



We then found the impulsive force that is consistent with 
the subsequent evolution of the system. Substituting the 
above results in the last equation in ([9]) we obtain a quite 
simple relation between the initial angle 0o and the ini- 
tial velocity 0o (where we suppress the argument 0+). 
Given the initial angle of the static disk, the angular ve- 
locity it acquires immediately after the driving apparatus 
is turned on is 



■ mdluj 

Po — — - — cos 00 



(10) 



Only pairs (0o,0o) related through this expression are 
valid initial conditions. We are now in position to show 
the relevant regions and curves in the space of initial 
conditions. In Fig. [21(a) we show the normal (blank) and 
reverse (dotted) regions in the space (0o,0o), separated 
by the synchronization line (grey curve). The black line 
represents the possible initial conditions as given by PU)) . 
We used the following parameters: uj = 6 rad/s (i/ ~ 1 
Hz), TO = 100 g, i? = 10 cm, I = 7 cm, and d = 10 cm. 
For these values, the constant of the motion that satisfies 
Eq. (O is f ~ 0.0543 J. We see that, for these parame- 
ters, no reverse rotations can occur (the black curve do 
not reach the dotted region). In contrast, if we set d = 40 
cm, keeping the other parameters (leading to £ w 0.2017 
J), we obtain the result displayed in Fig. [2Kb), where 
it is clear (the grey and black curves cross each other) 
that reverse rotations take place for an interval of initial 
angles centered in 0o = tt. The three initial conditions 
marked with a square, a circle, and a triangle represent 
the possible regimes: normal rotation, synchronization. 





FIG. 2: The dotted regions represent initial conditions leading 
to reverse rotations. We set d = 10 cm in (a) and d = 40 
cm in (b). The grey and black lines are the synchronization 
and initial condition curves, respectively. In (a) no reverse 
rotations can develop since the two lines do not cross. In (b) 
reverse behavior is possible for an interval centered in 0o = tt. 
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FIG. 3: Evolution of the initial conditions indicated in Fig. 2 
(b): square (black), circle (grey), and triangle (light grey). 



and reverse rotation. The time evolution of for these 
three conditions is shown in Fig. [31 The corresponding 
initial angles in radians are: 0o « 2.20 (black), 0o ~ 2.44 
(grey), and 0o ~ 2.68 (light grey). Notice that in the 
black curve, e. g., we have normal and reverse instan- 
taneous motions depending on the instant we record the 
velocity 0. The results presented in Fig. [2l refer to the 
global behavior of 0. 

It is possible to establish in a more precise way which 
initial configurations lead to reverse rotations. First we 
note that the range of initial angles that lead to reverse 
behavior is bounded by the intersections of the synchro- 
nization line ^ and the initial condition curve (HI 
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These boundary angles are given by 

2,(6) Ip ( 2g 2mdlw'^ ^ 

'''' =MM^1^7p-^; — " j ■ 

Since < cos^ (p^ ^ < 1, in order to get reverse rotations 
we must have 

2£ 1 

2mdl + Ip < ^ < — (Ip+ mdlf , (12) 

LU^ Ip 

where we have a condition imposed on £ = 
£(m, d, I, Ip,Lu). Let us analyze the limiting cases: ^ — 
2mdl + Ip and ^ = j^{Ip + mdl)^ . Substituting these 
expressions in the transcendental equation ([5]) we ob- 



tain relations involving a single parameter, K 




f VoH~T and K (^^^j = f (ct + 2), respectively, with 

a — 2mdl/Ip. The first equation has only the solu- 
tion tJ = 0, which is physically trivial, implying that 
the left-hand side of the inequality (fT^ is always ful- 
filled. The second equation, besides cr — 0, presents the 
solution a ~ 2.523. Thus initial configurations obeying 
2mdl/Ip — a, separate regions where reverse rotation is 
possible from regions where only normal rotations can 
occur. This condition involves only the relative scales 
D = d/R € (0, oo) and L — l/R £ (0, 1) and, in the case 
of a disk, reads 

£) = 0.631 + 1.261 L . (13) 

The above result is universal, in the sense that it is valid 
for any m and uj, and is independent from the absolute 
scales R, I, and d. As indicated in Fig. [H initial geomet- 
rical configurations located below the curve always 
lead to normal behavior, while configurations above it 
may enable reverse rotations, depending on the initial 
angle 4>o (the precise values of 0o are given by Eq. (Ilip V 
We also note that there is a value of d below which no 
reverse rotation occur. It is given by Dmm = ct/-\/2, that 
is, drnin — 1.784i?. For L > 0.5, variations in this param- 
eter produce virtually no change in D, which becomes the 
only relevant parameter to define the possible regimes of 
the system (see the "plateau" in Fig. S]). 

We stress that the obtained results are very general 
since we used the particular form of Ip only to obtain 
Eq. (fT3|) . For an arbitrary RB we have Ip = mP + 
^roR? , where 7 is a number and 7?, is a characteristic 
scale. For the disk we have 7 = 1/2 and TZ — R, while 
for a rectangular plate of sides a and b we have 7 — 1/12 
and TZ = v^oM-^. The general form of Eq. p3|) is 

D^^{jL-'+L), (14) 



with D = d/TZ and L = l/TZ, for TZ ^ 0. We note that 
in the "degenerate" case of a point mass connected by 
a massless rod to the pivotal point, obtained from the 
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FIG. 4: Configurations located above the curve may develop 
reverse rotations. Below the curve only normal rotation is 
possible. 



disk by taking i? = 0, we get a = 2d/l and the reverse 
condition becomes a one-parameter relation d/l > 1.261. 

We investigated reverse rotations in the circularly- 
driven motion of a RB whose intrinsic angular degree of 
freedom was shown to evolve according to a combination 
of pendular and uniform motions. This enabled the com- 
plete determination of the initial configurations that lead 
to reverse behavior (Eq. (fT2| V In addition a scale free, 
purely geometrical relation, defining the regions were re- 
verse rotation is possible was obtained (Eq. (fT4|) ). The 
effects of friction, assumed to be negligible in this work, 
may play an interesting role. The friction generated in 
the small contact area connecting the thin rod and the 
RB gives an extra torque (but no net force). This may 
be the only friction in the problem if we assume that the 
RB is kept in the horizontal position solely by the rod. It 
may also be of interest to consider the system immersed 
in a viscous fluid. This would make our apparatus very 
similar to a bioreactor for tissue growth [SJ . It seems that 
the influence of the regime of intrinsic rotation of the tis- 
sue construct (our disk) is yet to be analyzed. These 
pending issues are presently under investigation. 
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